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SECTION  I 
INTRODUCTION 


The  principal  problem  of  interest  in  the  present  investi¬ 
gation  is  the  same  as  that  treated  in  [1],  i.e.,  the  stress 
analysis  of  a  composite  laminate  built  of  anisotropic  elastic 
layers  of  uniform  thickness  and  subjected  to  prescribed  tractions 
and/or  displacements  on  its  boundary  surfaces.  The  body  is 
bounded  by  a  cylindrical  edge  surface  and  upper  and  lower 
faces  that  are  parallel  to  the  interfacial  planes.  This  assumption 
is  made  only  for  convenience  in  writing  the  governing  equations. 
There  is  no  difficulty  in  extending  the  model  to  include 
laminates  of  variable  thickness. 

In  practical  applications,  numerous  layers  may  be  present 
(use  of  100  layers  in  aircraft  structures  is  not  unusual).  Contemp¬ 
orary  models  are  incapable  of  providing  precise  resolution  of  the 
local  stress  fields  in  the  vicinity  of  stress  raisers  under 
such  conditions.  Global  models,  which  follow  from  an  assumed, 
usually  elementary,  displacement  field,  lead  to  the  definition  of 
effective  (or  smeared)  laminate  moduli  and  are  not  sufficiently 
accurate  for  stress  field  computation  [1].  On  the  other  hand, 
local  models,  in  which  each  layer  is  represented  as  a  homo¬ 
geneous,  anisotropic  continuum,  become  intractable  as  the 
number  of  layers  becomes  even  moderately  large  -  in  some  methods 
as  few  as  four  layers  result  in  technical/economic  barriers  to 
accurate  stress  resolution.  In  this  work  we  blend  these  concepts 
into  a  self-consistent  model  which  can  define  detailed  response 
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functions  in  a  particular,  predetermined  region  of  interest  (local) , 
while  representing  the  remainder  of  the  domain  by  effective 
properties  (global).  Such  dual  representations  are  not  without 
precedent  in  solid  mechanics.  For  example,  Gurtin  [2]  discussed 
this  approach  with  reference  to  the  solution  of  crack-tip  stress 
field  problems.  Wang  and  Crossman  [3]  used  an  effective  modulus 
representation  of  regions  of  a  laminate,  however,  only  the 
extensional  response  of  the  regions  were  considered,  i.e. ,  the 
flexural  and  flexural-extension  coupling  characteristics  of 
laminated  bodies  were  ignored.  Hence,  that  approach  fails  to 
provide  correct  solutions  to  certain  elementary  laminate  problems 
for  which  exact  solutions  are  available.  Stanton  et.  al.  [4] 
used  a  global  representation  based  upon  a  three  dimensional 
laminate  model  developed  by  Pagano  [5]  which  is  based  upon  the 
assumption  that  the  stress  field  is  only  a  function  of  one  space 
coordinate.  This  is  a  generalization  and  improvement  of  the 
material  model  given  in  [3],  however,  this  approach  is  not 
convenient  for  coupling  with  the  model  presented  earlier  [1]. 
Furthermore,  it  is  desirable  to  retain  the  model  [6]  as  a  special 
degenerate  case  of  a  global  model  since  that  result  was  shown  to 
produce  very  good  agreement  with  a  known  elasticity  solution  for 
transverse  normal  stress  a  [7]. 

Z 

There  have  been  several  investigations  of  the  interlaminar 
stress  fields  in  laminated  composites.  Pagano  [1]  has  given  a 
detailed  description  of  the  relevant  literature  in  this  field. 

A  recent  review  article  [8]  by  Solomon  presents  an  up  to  date 
literature  survey  in  related  topics  as  of  1980.  In  the  present 


2 


paper,  reference  will  be  made  of  only  those  publications  which 
are  not  covered  by  [1]  and  [8].  Spilker  and  Ting  [9]  have  con¬ 
ducted  the  static  and  dynamic  analysis  of  composite  laminates 
using  hybrid  stress  finite  elements.  Raju  et.  al.  [10]  have 
investigated  the  free  edge  stresses  in  layered  plates  using 
eight  node  isoparametric  elements.  In  both  these  publications, 
the  laminate  idealization  for  a  reasonably  accurate  finite 
element  analysis  had  to  be  very  fine,  i.e.  a  quarter  of  the 
laminate  was  divided  into  about  600  elements.  No  more  than  four 
layers  were  considered  for  numerical  calculations.  For  moder¬ 
ately  large  number  of  plies  (say  10) ,  these  approaches  will 
lead  to  computer  storage/economic  difficulties. 

Blumberg  et.  al.  [11]  studied  the  edge  effects  and  stress 
concentrations  in  composite  laminates  made  of  glass  sheets  bonded 
with  polymer  adhesive.  The  governing  equations  employed  were 
similar  in  nature  to  those  given  in  [1] ,  however,  not  as  general. 
For  example,  only  isotropic  layers  were  considered  with  the  stiff 
layers  being  represented  by  the  Kirchhof f-Love  theory.  Further¬ 
more,  the  implied  edge  boundary  conditions  are  not  sufficient  to 
satisfy  the  principle  of  "layer  equilibrium"  [1] .  The  differential 
equations  were  solved  by  perturbation  technique  defining  the 
dependent  variables  at  three  different  regions  along  the  width  of 
the  laminate.  This  division  of  the  width  has  enabled  the  authors 
to  overcome  computation  overflow/underflow  difficulties. 

Finally,  Partveskii  [12]  has  presented  an  approximate  treat¬ 
ment  of  a  free  edge  problem.  This  model  combines  the  treatment  of 
[13]  with  a  model  based  upon  a  layer  on  an  elast^  found  .ion  in 

order  to  define  the  distribution  of  interlaminar  normal  stress. 
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SECTION  II 

VARIATIONAL  PRINCIPLE 


The  laminate  considered  in  the  present  investigation  is 
shown  in  Figure  1(b).  The  laminate  thickness  comprising  of  N  + 

M  layers  is  divided  into  two  parts  viz;  (i)  local  region  (£ ) 
and  (ii)  global  region  (g) .  N  is  the  number  of  layers  in  the 
local  region  and  M  is  the  number  of  layers  in  the  global  region. 
In  this  work,  we  shall  assume  that  the  interface  between  g  and  JL 
is  a  plane  2  =  constant,  although  less  restrictive  assumptions 
are  possible.  A  variational  principle  as  described  below  has 
been  used  to  derive  the  set  of  field  equations  and  boundary 
conditions.  Different  variational  functionals  in  two  different 
regions  of  the  laminate  are  used  such  that 


and  body  forces  are  neglected. 

In  equation  (1),  the  first  term  represents  the  potential  energy 
for  the  global  region,  the  second  term  is  the  Reissner  varia¬ 
tional  functional  for  layers  in  the  local  region  and  the  third 
term  is  the  potential  energy  of  the  prescribed  surface  tractions. 
The  notation  used  here  is  the  same  as  that  of  reference  [1],  with 
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the  exceptic"  that  the  subscripts  ' l '  and  '  g’  denote  respectively 
the  local  and  global  regions.  In  equations  (2)  ,  w  and  w  are  strain 
energy  density  functions,  the  first  in  terms  of  displacements  u^ 
and  j »  the  expansional  strain,  components,  and  the  second 
in  terms  of  stresses  and  e^^  j  . 

For  a  layered  continuum  in  the  local  region,  equation  (1) 
can  be  rewritten  as 


r  N 

/  [5  *  uj,i*  •  "lll°dvk  -  J,  "ids  + 

U-i  v,.  s 


wdv  >  =  0 . 


where  the  superscript  (k)  attached  to  the  bracket  signifies 
that  each  variable  within  the  bracket  is  associated  with  the 
kth  layer.  The  use  of  Green-Gauss  theorem  and  some  mathematical 
manipulations  in  equation  (3)  yield  the  following  equation. 


6a.  .  -  a.  .  .  5u.  ]  (k)dv. 
3  cri  j  13  13  0  i  k 


E  /  ‘Oij  -  oij,j«ui),k)dvk 

k=1  v 

r  r  n-1 

J  (t.-  x.)6u.ds  +  J  t  .  6u.  ds  +  V  f  (xfk)5uf 
s'  1  s"  Z-J  1  1 


k=i  I, 


W'tl)>dik-  / 

vg 

Tft)+xfg)  )  6u .  ds  =  0 

X  11 


a . .  . Su . dv  + 

ij  '3  i 


/ 


( T^-T^ ) 6uids  +  J  Ii£uids 

s" 

g 

(4) 


Vk  is  the  volume  of  kth  layer,  and  S£  represent  the 
■  surfaces  bounding' the  local  region,  the  former  representing 
>ortion  with  prescribed  tractions  and  the  latter  with 

II  I 

:ribed  displacements.  1^  represents  the  interlaminar 

ice  between  kth  and  (k+i)th  layers  in  local  region  that  does  not 
/{ 

tg  to  S'  or  S  ,  S'^  and  represent  the  bounding  surfaces 

le  global  domain  with  prescribed  traction  and  prescribed 
icement  conditions,  respectively.  S  represents  the  surface 
i  to  the  local  domain  and  the  global  domain.  A  superscript/ 
ript  l  denotes  the  local  region  and  g  denotes  the  global 

Clearly,  as  shown  by  equation  (4) ,  the  governing  equations 
asticity  can  be  obtained  as  a  consequence  of  variational 
Lon  (1).  Equation  (4)  will  be  used  to  derive  the  field 
ions  and  boundary  conditions  for  the  two  regions  of  the  laminate. 
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SECTION  III 


DEVELOPMENT  OF  THEORY 

For  each  layer  in  the  local  domain,  the  theory  developed  in 
[1]  has  been  used.  The  details  of  the  derivation  of  equilibrium 
equations  and  continuity  and  boundary  conditions  for  this  domain 
are  not  repeated  in  this  paper.  For  the  sake  of  continuity,  only 
relevant  equations  are  provided.  Figure  2  shows  the  coordinate 
axes  and  thickness  of  a  single  layer  in  the  laminate.  The  inter¬ 
laminar  stresses  cr,  r  „  and  t  at  the  top  of  the  layer  are 
denoted  by  p2,  t2  and  s2,  respectively,  while  the  corresponding 

stresses  at  the  bottom  of  the  layer  are  designated  as  p-^,  t.  ,  and 
s^.  In  the  local  domain  the  inplane  stress  components  are  assumed 
to  vary  linearly  through  the  thickness  of  each  ply.  The 
substitution  of  these  stress  components  in  the  differential 
equations  of  equilibrium  [1]  yields  the  interlaminar  stress 
components  in  terms  of  tractions  p.^,  t^,  si  (i  =  1,2)  and  force 
and  moment  resultants.  These  stress-stress  resultant  relations 
have  been  used  in  equation  (4)  in  local  domain  integrals.  In 
the  global  domain,  an  assumed  continuous  thickness  distribution 
of  the  displacement  field  is  used.  On  the  basis  of  these 
facts,  the  field  equations,  interfacial  boundary  conditions, 
and  edge  conditions  within  the  local  continuum  remain  the 
same  as  derived  in  reference  [1].  The  development  of  the 
required  relations  for  the  global  domain  and  global/local 
interface  follows.  We  assume  that  the  global  domain  is 
composed  of  layers,  each  possessing  a  single  plane  of  elastic 
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symmetry,  z  =  constant .  In  this  domain,  the  displacements  are 
assumed  to  be  of  the  form 


u  =  u° (x,y)  +  zY  (x,y) 

v  =  v°(x,y)  +  zV^(x,y)  2  ^ 

w  =  w°  (x, y)  +  z'i  (x,y)  +  o(x,y) 

where  u,  v  and  w  are  the  displacement  components  in  the  x,y  and 
z  directions,  respectively.  It  can  be  seen  that  the  number  of 
displacement  functions  agrees  with  that  given  by  the  variational 
principle  for  each  layer  in  the  local  domain.  The  substitution 
of  the  displacement  functions  (5)  into  the  strain  displacement 
relations  of  elasticity  leads  to  the  following  stress-strain 
relations 

°i  =  Cij  (ej  +  z*j  '  ej}  <*0  =  1,2, 3,6) 

°i  =  Ci3  (Ej  +  Z*j  +  T  V  (i'3  "  4*5) 


in  standard  contracted  notation,  where  are  the  components  of 
the  anisotropic  stiffness  matrix,  ej  are  engineering  expansional 

O 

strain  components  and  c^,  Ki  and  8i  are  defined  by 


0  _  3u 

ei  ax  ' 


ro  _  9v° 
£2  ay 


=  p 


e°.  =  V  +  3w 


y 


By-'  "5 

3* 


°  ,  .  9w 

£_  =  Ip  + 


Z 

O 


ax  ' 


°  _  3u  av 


ay 


ax 


*1  axX'i-  *2  ~  Jy~'  *3  ~  * 
bpz  ftp-  a p  „ 

*4  =  dy~'  *5  =  TT’  *6  =  W  +  33? 


(7) 


ft  -  JL£  r  _  34> 
e4  “  3yr  S5  "  ax 


The  stress  components  a^,  a 3,  <?4, 

°2 ,  Tyz,  txz,  txy,  respectively  while 


ac,  a.  stand  for  a  ,  a  , 

5  6  x  y 

O 

(i  -  1,2, ...6)  represent 


p  1 


corresponding  engineering  inplane  strain  components.  We 


introduce  the  following  stress  and  moment  resultants 

H/2 

N ^  -  J  a^dz 
-H/2 
H/2 

=  y*  zdz 

-H/2 


(i  *  1,2, 3,6) 
(8) 


where  H  is  the  thickness  of  the  global  region  and  N3  and  m3 
are  mathematical,  not  physical,  quantities. 

Substituting  equations  (6)  into  (8)  and  conducting  the 
integration,  we  obtain  the  following  constitutive  relations 

for  the  global  domain 


Na  AaBe6  +  Ba8*B~Ra 


M  =  B  _e®  +  D  0k  -Q 
a  aB  B  aB  8  a 


(ot,e  =  1,2, 3,6) 


VA.  WW  1 

V.  =  A. . e?  +  B. +  4  D. .8. 
i  13  3  13  3  2  13*3 


R.  =  B.  .e?  +  D.  .* .  +  ^  P.  .8  . 
i  13  3  13  3  2  13  3 


(i,j  =  4,5) 


where 


s.  -  i  D. .e?  +  i  F. .k .  +  t  H..B. 
1  2  13  3  2  13  3  4  13  3 


(A*,  .,'b.  .,  D,.,  F..,  H .  . )  =  /  (l,z,z2,z3,z4)  C.  .dz 

13  13'  13'  13'  13  J  n 


-H/2 


2  / 1  •  •  ,6) 


-H/2 


C«8e8dZ 


Qa  = 


/  2C«6eedJ 

-H/2 


(a, 8  *  1,2, 3,6) 


(Vi,Ri,Si)  *  rH/2ai(l,z,z2)dz  (i*4,5) 


/ 

-H/2 
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With  the  knowledge  of  the  distribution  of  elastic  properties 
and  expansional  strains,  one  can  obtain  the  values  of 

^  A—' 

effective  stiffness  matrices  A,  B,  D,  P  and  H;  and  effective 

"nonmechanical"  stress  and  moment  resultants,  P.  and  5.. 

1  1 

As  in  [1],  we  make  the  following  definitions 


/ 


H/2 


-H/2 


f  ( 1  — 
ItA,  H' 


2dz 

H 


(10) 


where  f  may  represent  any  of  the  displacement  variables  u,  v 
and  w.  Through  the  use  of  relations  (5)  and  (10),  the  functions 
involved  in  the  displacements  for  the  global  domain  can  be 
expressed  as 


u  = 


r- 


V  = 


K  = 


I- 


*  1  w* 

*  -  4  <S  - 1) 

H  J 


(ID 


With  these  relationships,  one  can  express  the  constitutive 
relations  for  the  global  domain  in  terms  of  the  same  displace¬ 
ment  parameters  as  those  for  the  local  domain.  This  simplifies 
the  definition  of  the  required  continuity  conditions. 

With  the  use  of  the  assumed  stress  field  in  the  local  region 
and  displacement  fields  in  the  global  region  in  equation  (4), 
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the  required  field  equations,  continuity  and  boundary  conditions 
can  be  obtained.  For  the  local  region  the  equilibrium  equations, 
constitutive  relations,  edge  conditions  and  interfacial  continuity 
conditions  are  given  in  II].  For  the  global  region  the  equili¬ 
brium  equations,  which  follow  from  substituting  (5)  into  (4) , 
become 


where  the  symbols  t^,  s^,  and  (i  =  1,2)  retain  the  same 
meaning  as  defined  earlier  for  the  local  domain.  Assuming 
perfect  continuity  of  tractions  and  displacements  at  the  g-Jt 
interface,  the  local-global  interfacial  conditions  are  given 
by  the  previous  substitution  into  (4),  as 


(k)  _  . (k+1) 
2  ~  rl 


(k+1) 

S1 


(13) 


11 


r 


and  (a)  local-global  (interface  I  of  figure  la) 


1S5-S45VS55T51<k)  “  l-“*  *  f  *xI(k+1> 


<k> 


=  [-v°  +  #  *.] 


4  444  45  5 


^y2“S33R2^  (  5  -  r“ 


(k+1) 


[-w‘  +  |*2  -^^](k+i) 


2  ryJ 
2 
8~ 


(13a) 


It  can  be  shown  that  if  we  consider  more  than  one  global  region, 
the  following  interfacial  conditions  are  required: 

(b)  global-local  (interface  II  of  figure  la) 


[u"+|*x](k)  .  t“3-s45Q4-s55Q5] (k+1) 


[v°+fv 


(JO  _ 


"  [ot4-S44Q4"S45Q55 


(k+1) 


(13b) 


o.H„.  j_H  (k)  _  r„  „  t  (k+1) 


[w  WV1 


=  [y1-s33r1i 


(c)  global-global  (interface  III  of  figure  la) 


o » M,i.  ,  (k)  _  i  (k+1) 


[U0^1 


-  [n°-rx] 


(13c) 


(v°+§V  (k)  =  [V°"f^y3  (k+1) 

[w°+^z+S  <j>]  (k)  =  cj> ]  (k+1) 

where  the  parameters  with  superscript  'k+1'  represent  those  for 
the  layer  above  the  kth  layer.  The  parameters  on  the  left-hand 
side  of  equations  (13a)  and  right-hand  side  of  equations  (13b) 
are  defined  in  reference  [1],  In  the  expression  for  R2  of  [1] , 
the  roles  of  and  p2  were  interchanged  by  mistake  whereas  R^ 
was  correct.  The  correct  expression  for  R2  is 


R2  - 


(6p2+p1)h/-7hNz-30M! 

73h  f 


(14) 


For  the  edge  surface  of  the  global  domain,  one  term  each  of  the 
following  products  must  be  prescribed: 
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W  Nnsus'  Vn'  MnsV  V  '  Vz'  V* 


(15) 


The  boundary  conditions  at  the  top  surface  are  given  by 


(N+l)  _  ?(N+1) 


or  u°  +  |  "  V 


(16) 


(N+l)  „  -<N+1)  or  vo  +  H  ^  =  ~ 


*2 


(N+l)  .  ~(N+1)  or  ^  +  H  ^  +  £  * 


W 


where  the  right  hand  sides  in  the  aforementioned  equations  (16) 
represent  the  prescribed  external  tractions  or  displacements. 

The  boundary  conditions  at  the  bottom  surface  remain  the  same 
as  those  explained  in  reference  [1].  This  completes  the 
development  of  the  present  theory.  We  observe  that  the  governing 
equations  for  the  global  continuum,  equations  (9),  (12-13), 
combined  with  the  governing  equations  for  the  local  continuum, 
equations  (25-28)  of  [1],  and  boundary  conditions  at  the  bottom 
and  the  top  surfaces  constitute  a  set  of  23N  +  27  equations  in 
terms  of  like  number  of  unknowns.  This  system  can  be  reduced 
to  13 (N+l)  equations  by  eliminating  the  force  and  moment 
resultants  from  the  set  of  governing  equations.  Relations  (15) 
show  that  7  edge  conditions  are  required  for  the  global  domain, 
while  7N  edge  conditions  are  required  for  the  local  domain, 
equations  (29)  of  [1J. 
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SECTION  IV 
PROBLEM  DESCRIPTION 

The  present  model  has  been  used  to  conduct  the  free  edge 
stress  analysis  in  a  symmetric  laminate  consisting  of  2 (N+M)  per¬ 
fectly  bonded  layers,  see  Figure  1.  The  laminate  is  subjected  to 
forces  applied  only  at  the  ends  y  =  const,  such  that  a  constant 
axial  strain  =  e  is  imposed.  Because  of  the  symmetry  in  ply 
orientation  of  the  laminate  about  the  midplane,  the  deformation 
is  symmetric  with  respect  to  x  and  z.  Only  the  z  symmetry  will 
be  employed  in  the  specific  problem  treated,  so  that  half  of  the 
laminate  thickness  has  been  considered.  The  lower  N  layers  form 
the  local  region  whereas  the  remaining  M  layers  constitute  the 
global  region.  The  stress  field  in  this  class  of  problems  is  a 
function  of  x  and  z  alone  and  consequently  the  force  and  moment 
resultants  and  the  interlaminar  stresses  depend  on  x  only. 
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SECTION  V 
SOLUTION 

1 .  LOCAL  DOMAIN 

By  the  use  of  strain  displacement  relations  (1)  of  [14] 

it  can  be  shown  that  the  most  general  form  of  the  weighted 

displacements  within  each  layer  is  given  by 

C.y2 

u  =  U  (x) - -  +  C3y 

v  =  V(x)  +  C^xy  +  C2y 

hw  =  W (x)  -  6C,.xy  +  3C4y  -  3Cgy2 

u*  =  ip  (x)  +  C5y  (17) 

v*  =  fi(x)  +  Cgy 

w*  =  <f(x) 

/\  2 
hw  =  x(x)  -  2C$xy  +  C4y  -  Cgy 

where  U,  V,...  x  are  arbitrary  functions  of  x  and  C.  ,  C_...C,. 

L  c  o 

are  constants.  We  should  recall  that  the  stress  strain  relations 
(3)  of  [14]  must  be  written  for  each  layer.  The  use  of  the 
foregoing  relations  in  the  strain-displacement  relations  yields 

el  =  I  U'x  =  I0’  (x) 
e2  =  7  v,y  =  |(C1x+C2) 

^  =  3  w*  =  3<^(x) 

e6  "  I  <5>y+"<x>  *  7  [V  (x)+C3l 
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I 

\ 


(18) 

<3  =  |  h  (3w-w)  =  |  [3X  (x) -W (x)  ) 

2  2 
<6  “  V  (u*'y+V*'x>  =  V  ^’(X>+C5] 

e4  =  x  <w/y-w,y)+^-  =  |  [2^(x)-2C5x+C4] 

es  =  T"  ("'x“'''x)+rr =  I  [w' (x)-x'  (x)+4*(x) ^ 


in  each  layer,  where  h  is  the  layer  thickness.  The 
substitution  of  the  strain  components  (18)  in  the  equilibrium 
equations  (28)  of  [1],  through  the  constitutive  relations 
(25),  [1],  gives 

AlllU',+3A13*'+A16?V''+A13Tr  (P1'x*P2'x)+t2-tl-Aljej.x-^rhCl 


A16lU"  +  3A36*'+A66^"+^^T^  (Pl,x+P2.x,+»2-sl=A6je1,x-Tihcl 


36  662 
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A135U,+3A33<’+A367V'1ro  <A33s33*5,lpl+P2,41ij(tl,*-t2.x1  * 


A  3 jej”A2  32 (C1X+C2  > ”A36JC 3 

Bllh2  S  B  h2  B13S33h2  5 

nr  *"+  !b13[3X‘-w,]+-1|-  a"-.  (P2>x-plfX)-fA45^ 


“fA55 (w’- X,+4^)+|§(t i+t2)  =  |A45(C4-2C5x) 

B1fih2  «;  BK,.h2  B.fiS..,h2  = 

-if-  ♦■♦§B36(3*'-»’H-T-  <P2.x-Pi,x)-K4[! 


16 


~^4  ^  (w‘  -  X*  +4  Jj)  +yy  (s^+S2>  -  |*44(C4-2C5X) 
§A45fl'+|A55(w--X"+4*')1|(tli)t+t2>x)+p2-Pi  =  |a45C5 

g 

Ks*'^  33S33-14'  <P2-Pl>-II(tl,x+t2,x>i 


"fB23C6'lB36C5 


(19) 


These  equations  are  valid  within  each  layer.  Similarly  the 
remaining  field  equations,  the  interface  continuity  conditions, 
are  given  by 


{  §X'-5«,-J*,-U-|*-||[s55(3t1-t2)+S4s(3Sl-S2)l]  (k+1) 

+  (l’<'-IW'-4>,-U-|*,-A1S55(3t2-tl)+S45(3s2-sl)1}  <k) 

+C5  (k)  )y+(C3(k)-C3(k‘*'1)  )y+|(C1  (k+1)-C1  (k)  )y2. 


=  3  (r  CK) 


V“lQ'Il[S45(3tl“t2)+S44(3srS2)1l  (k+1> 


(■ 

✓ 

+  ^-V-|fi-^[S45(3t2-t1)+S44  (3s2-  Sx) 


(k) 


(k+l>  „  (k) 


=  f(C5^'i'+C5'7V/)x-|(C4'Jvy+C 


(k)  ,  (k+1) 


) +3 (C. 


{k+1)+C#.  (k>  )y 


+  (C1(k)-C1(k+1) )xy+(C2 (k)-C2 {k+1) )y . 


{•% 


35 
225B 


7A-  «h  i  a  p  7A  - -  h  ^  p 

-f1  U,  +  (-2lA33+^f)(t  -  -r-  V,+TB13h^ 


+  (- 


_11  _  525  )  +  ,  25  +  105  .  15B36h  fl1 
2h  2S^hJ  X  +  (  2hB33+2S33h)W^  2  0 


15B33S33 


7A 


15B 


+  hPl<6-7-ir  S33-  14  JJ)+hP2(1-lF  S33+-lfi  S33> 


) 


(k+1) 
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225B,,  cot;  -7C  me  15B--h 

+  2^h)x  +  (iB33-^3h)W— 2^  «* 

7A33B33  ^B33B33  ^A33B33  B33 

+hpi(1 — n-^  +  s33 


(k+1)  (k) 


(k)  (k+1)  (k+1)  (k) 


^1  (k) 


=  6<i 


hf  )xy+3(hf  - 


C.  C  ,  o 

4  \  . o  /  6  6.2 

r —  )y+3(r - ir-)  y 

hk+l  hk+l  hk 


(k+1)  (k+1)  (k+1) (k) (k)  (k)  , 

/  r>  it  /-t  .  r«  «  L  A  v  .  A 


"10 (S33A23hk+lCl  +S33A23 


(k)  i  r  r 

hkCl)+^(S33/ 


^9  ^ 

7- r -  V’-r-  C3+7A3jej 


■TB23hC6  + 


15B  h  15Boch 

^  J  p  JO  r* 

5  6  2  5 


1_5_  hc  l^(k+1)+-L  fs  f  7^11h r  7A36i 

2  B36hC5jJ  35  [33  “7  2  hC2-7  2 


■hC_+7A-  .  e  .  - 
3  3]  ] 


_  4.  (k+1) 


_  „  (k+1) 


(k)  _  (k+1) 


The  last  six  equations  are  valid  for  k  =  1,  2,  —  ,N-1,  since  we 
recall  that  N  is  the  number  of  layers  in  the  local  region. 

Due  to  the  symmetric  lamination  geometry,  the  interlaminar 
shear  stress  components  and  the  transverse  displacement 
component  W  all  vanish  on  the  central  plane,  z  =  0.  We  shall 
take  advantage  of  these  conditions  by  considering  only  the 
upper  half  of  the  laminate,  i.e.  z>0.  Thus,  the  boundary 


conditions  at  the  lower  surface  reduce  to 
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(1) 

3^  l— ^hU'  +  (-21A33+|^)<H^^hV'-H^B13h4i' 


^a'+hPl(6-7^jfc 


15B33S33}  (1-7A33  s  +15B33£  } , (1) 

14  +np2 11  1Q  S33+^4  ®33 


t6C_  3 C.  3 C,  , 

--h^xy+-H4  y-nry' 

+  TB23hC«+¥B3«hC5> 


2  x  B33  "^^2  3  "^A-- 

+T0  S33A2  3hCl"35~C  2  hC2  r^C3+7A3jej 


t^1J  =  s^1*  =  0  (21) 

Since  equations  (19-21)  must  be  satisfied  for  all  values  of 
y,  it  follows  that 

c4k)  =  C£k)  =  C^k)  =  0,  k  =  1,2, ...N 

and 

_(k+l)  _  _(k)  r(k+l)  _  r(k)  r(k+l)  _  _(k)  _  -  „  , 

so  that  C^,  C2  and  C3  are  the  same  for  each  layer.  Further, 
it  has  been  shown  in  [1]  that 

C1  =  c3  =  0 

C2  =  2  e  . 

We  will  use  these  values  of  in  the  forthcoming  work. 

2 .  GLOBAL  DOMAIN 

With  the  aid  of  relationships  (11)  ,  we  can  obtain  the 


following  strain  displacement  relations  for  the  global  domain 
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^  45  .a  w. 

*3  *  ^  (u  '  3> 


k6  ‘  I  (u*'y  +  v*'x> 

-  A 

re  u,  u*,  -  w  now  refer  to  the  global  domain.  We  may 

d  observe  that  the  relations  (17)  are  also  valid  for  the 
bal  domain.  The  substitution  of  the  values  of  inplane 
ains  and  curvatures  in  the  stress-strain  relations  (9) 
ough  equations  (17)  and  their  subsequent  use  in  equations 
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(12)  yield  the  following  form  of  equilibrium  equations 


I^uu"  +  1*13  ♦'  +  5*16  V  +  I V  *  ih. 


.  45^ 

+  "3  B13 
H 


<X'  -  1  w'>  +  *2  -  *1  ‘  Pl.x  *  0 


2  61 


IX,  o-  +  |  a63  +  |\.fi  v-  +  1  b61  r  +  -  b66  n" 


♦  ^^63  <*’  “  |  W-)  +  s2  -  sx  -  P6  x  =  0 


1 

1  "31 


A„  0*  -  |  D55  4”  *^33**!^  •'’?  ''31  ”  "55' 


A/  Aj 


V  +  ~  (B*.  -  B--)  i> 


♦  \  ‘*36  *  <J1  *55  -  %  F55>  X"  *  ^33  X 


(i%  -  4r  B55)  W»  -  ^  B33  W  ♦  A32  E  -  P3  -  5  *Pl 


15  r~  —  H 


-  T  (Pi  +  P7)  =  0 


2  H 


3  55 


H 


1-  3*. .  U"  +  t,  +  I  B16  V”  +  H  DH4>”  "  H  A55  ^ 


2  I311U  H  '  13  55 


+  I  °16  -  1 *45  fl  +  (7  °13  +  F|A55  "  ^3  °55>  *’ 


45 


,15  n  .  JL  a 
"  (_3  D13  8H  55 
H 


15D 


2h" 


|^)W'  +  f  (t2  +  tL)  -  Q1#x  -  0 


i  B„  U"  ♦  |  (¥63  -  B^)  *  T  <>66  V"  +  5  °61  -  if  A45  * 


2  61 


+  1  °66  *  I  A44  11  +  (T  D63  +  S§  A45 


45D 


45 


)  X’ 


2H" 
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N|H  +  tfi|U> 


0 


-  3  °63  *  ¥h^45  -  3  D45>  W'  +  !  (S1  +  s2»  '  °6,x 

H  2H 


B55  4>"  +  H  A55  +  H  A45  fl’  +  ^^3  °55  ’  8H  A55)X” 


,  9  15 

8H  A55  '  2fJ3  "55 


D«)  W"  +  p  -  p.  =  0 


Bn  -  -  f55 


>"  +  \  ^33  ♦  +  1  B36  v‘  +  I  (D3r^i)  *'  (23) 


H  '  31  2 


D, 


+  \  **-3'  C'  +  <Mh  °55-^  H55>  X"  +  p  °33  « 


+  <7^  H55-tIh  D55>  w"  -  3  °33  ™  +  B32  £  *  «3  '  V  “W  *  0 
4  H  H 

As  for  the  local  domain,  in  the  foregoing  equations  we  have 
used 

C1  =  c3  =  c4  =  C5  =  c6  =  ° 


C2  =  2e 


In  equation  (23)  the  effective  modulus  matrices  defined  in 

equations  (9)  are  used.  The  continuity  conditions  (13) and  (13a) 

at  local-global  interface,  on  substitutions  of  (17)  and  (14)  of 

the  present  paper  and  equations  (16)  of  [1],  reduce  to 

. (N)  _  (N+l) 

c2  ~ 


(N)  _  (N+l) 

s2  S1 


(N)  _  (N+l ) 
P2  "  P1 
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(f  x'  ■ ■>'  +§♦■ '  -U-f*-l2  Is55  ( 3t2'tl)+S45<3s2'sl)  j}‘N> 


=  (-U  +  3  4i) 


(N+l ) 


{-H* 


X2^S45(3t2"tl)  +  S44(3s2 


,  1‘\(N) 

-Sin] 


=  (-v  +  3  n) 


(N+l) 


{  15i[-^h0 1  +  <-21A33+S^)  «-IA36hV ’  -¥hBl  3* ' 


+<-  *  (5KB33-5if;»M-irlB36i1’ 


+  ho  (1  7A33c  | 15B33S33>  h  (6  7A33S33  15B33S33x , 

+  hp1(l-  1Q  S33+  14  +hP2  6  1Q  n  )] 


(N) 


M  x  - 3  *  -  Aw 


(N+l) 


(25) 


The  boundary  conditions  at  the  top  surface  considered  in  the 
present  investigation  are: 


(N+l)  _  . (N+l)  _  (N+l)  _  . 
s2  ~  t2  ~  P2  0 


(26) 


3.  EDGE  BOUNDARY  CONDITIONS 

We  now  turn  our  attention  to  the  edge  boundary  conditions, 
which  require  consideration  of  N  ,  N  ,  V  ,  M  , M  ,  t,  and  t., 
for  each  layer  on  x  =  +  b,  since  no  displacement  edge 
conditions  are  involved  in  the  present  class  of  boundary 
value  problems.  However,  all  these  functions  cannot  be 
independently  prescribed  because  of  the  consequences  of  interface 
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continuity  and  overall  equilibrium  of  the  entire  laminate. 


That  is,  the  interface  continuity  conditions  given  by  the 

00 

fourth  of  (20)  prohibit  arbitrarily  prescribed  values  of  t. 

(k)  (1)  (N+l)  1 

and  t2 .  Furthermore,  t^  and  t2  have  already  been 

specified  by  (21)  and  (26)  for  all  values  of  x.  These  relations, 

in  conjunction  with  the  first  equilibrium  equation,  see  (26) 

of  (I] ,and  equation  (12), can  be  used  to  establish  the  result 


N+l  (k) 

y  n  =  o 

X,  X 

k=l 


which  requires  that 

N+l  (k)  N+l  (k) 

Z  Nx<b>-  Z  Nx(_b)  =  0  (28) 

k=l  k=l 

(k) 

Therefore,  only  2N+1  values  of  N^  can  be  arbitrarily  prescribed 

on  the  edges  x  =  +  b.  We  can  make  the  same  statement  regarding 
(k) 

N  since  an  equation  of  the  form  (28)  can  be  derived  in 
xy 

similar  fashion  for  this  function.  Hence,  the  edge  boundary 

conditions  for  the  local  domain  may  be  expressed  as 

(k)  (k)  (k)  (k)  (k) 

Nx(+b)  =  Nxy  (+b)  =  Vx(+b)  =  Mx(+b)  =  (+b)  =  0 


(k) 

t2  (+b)  =0  (k  =  1,2,...  N)  (29) 


while  those  for  the  global  domain  are 


(N+l)  (N+l)  (N+l)  (N+l) 

Nx  (b)  =  Nxy  (b)  =  vx  (±b>  *  Mx  (±b) 
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(N+l)  (N+l )  (N+l ) 

-  M  (+b)  =  R  (+b)  =  Sv  (+b)  =0  (30) 

xy  x  x 


Thus,  the  present  boundary  value  problem  consists  of  the 
differential  equations  (19,  20,  23,  and  25)  subject  to 
boundary  conditions  (21,  26,  29  and  30). 


The  general  solution  for  each  dependent  variable  consists 
of  the  sum  of  two  parts:  i)  a  complementary  solution  defined 
by  the  homogeneous  form  of  (19,  20,  23  and  25),  and  ii)  a 
particular  solution.  In  the  particular  solution  (denoted  by 
subscript  P) ,  the  only  nonvanishing  functions  are  given  by 


(k)  (k)  (k)  (k)  (k)  (k) 

<)>  =  a1  Xp  =  +  a<>  Wp  =  ^a2  '  k=l ,  2 , . . .  ,N+1 

p  (3: 

(k) 

where  ai  (i=l,2)  are  constants  given  by  substituting  (31)  into 
(15)  and  (20)  to  get 


(k)  1  (k)  (k)  (k) 

al  =  TTTk)  (A3gVA23hke)  '  a»  =  °'  k=l,2,...,N 

ja33 

(l)  (1) 
a2  =  hlal 


(k+1)  (k)  (32) 

a_  (k+1)  (k)  a- 

; -  —  =  a,  +  a.  +  7-=-  ,  k  =  1,2...  N-l 

hk+l  1  1  hk 


where  aQ  =  he  (6  =  1,2, 3, 6). 


Similarly  using  (31)  in  (23)  and  (25),  we  get 


(N+l) 

al 


1 

~  (N+l) 
A-,-, 


3 


CA23) 


(N+l) 


an  = 


45 


^33^3~eA23^  ”  A33^3  “  eB23^ 


/  (b33-a33d33) 


(N+l) 


(N+l)  _  5 

a  2  "2 


ao 


+  H 


(N) 

a»  (N+l) 


(33) 


Since  the  field  equations  are  linear  differential  equations 
with  constant  coefficients,  the  complementary  solution  (sub¬ 
script  H)  for  each  dependent  variable  consists  of  a  series 
of  terms  of  the  form 

(k)  (kK 

fH  *  F  e  k  =  1,2 - N+l  (34) 

(k) 

where  f  represents  any  of  the  dependent  variables  and 
F^)  are  constants.  The  substitution  of  (34)  in  the  homo¬ 
geneous  form  of  equations  (19,20,  23  and  25)  yields  a  system 
of  13 (N+l)  linear  algebraic  equations.  This  set  of  equations 
can  be  written  in  the  following  matrix  form 

[J]  [P]  =  0  (34a) 

where  J  is  a  13 (N+l)  x  13 (N+l)  matrix  dependent  upon  material 

properties  of  the  laminate  and  X,  and  F  is  a  vector  of  constants 
(k) 

F  defined  in  (34)  .  For  a  nontrivial  solution  of  equation  (34a)  , 
the  determinant  of  the  coefficient  matrix  J  has  been  equated  to 
zero  and  the  resulting  equation  has  been  solved  for  specific 
values  of  X  by  the  method  of  Jenkins  and  Traub  [15] . 
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The  algebraic  expressions  for  the  elements  of  J  were  not  written 
owing  to  the  complexity  of  such  expressions,  even  in  the  simplest 
cases.  The  computer-calculated  values  of  X  have  been  used  to 
carry  out  the  required  analysis.  The  extreme  (highest  and  the 
lowest)  powers  of  X  in  the  polynomial  expansion  of  [J]  were 
investigated  for  small  values  of  N+l  and  deduced  for  arbitrary 
values  of  N  and  M.  It  was  found  during  numerical  calculations 
that  for  incompatible  extreme  powers  of  X,  the  solutions  for  X  were 
nonconvergent  in  the  iterative  procedure,  whereas  for  compatible 
powers  these  converged  very  rapidly.  Following  the  procedure 
described  in  [1],  the  following  observations  were  made  i)  only 
even  powers  of  X  occur  in  the  determinant,  ii)  the  lowest 
power  of  X  is  X4,  and  iii)  the  highest  power  of  X  is  x12^N+1^. 

As  in  [1],  the  functions  corresponding  to  the  repeated  zero 
roots  for  X  can  be  written  in  the  following  form: 

U„k)  =  A,  x  +  A 
Ho  1  O 

V^)  =c  x  +  C 
HO  1  O 


A(k)  _  _(k) 
♦ho  -  Bo 


(35) 


x(k)  .  E(k) 
XHo  O 


=  3(E*k)  +  d  ) 

HO  O  O 


dQ  =  o  for  k  =  1,2,...,  N 
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where 


(k) 

B 


(k)  (k) 

(A13A1  +  hk 


6A 


(k) 

33 


,  k  =  1,2. . .  N 


(1) 

E 


(1) 
-  hlBo 


(36) 


(k+1)  (k) 

E  (k)  E  (k+1) 

=  B  +  —■  +  B 

hk+l  o  hk  o 


k  =  1,2. . .  N-l 


(N+l) 

Bo  -  6(S-33D33-^) 


-H 


(A13D33“B13B33)  A1  +  (A36D33_B36B33 )C1 


do  15D 


«!_  h  <N 

D,,  )B33Bo 


(N+l) 


33  'v. 


+  |  (B^  +  B^  C,) 


(N) 

(N+l)  (N)  E  (N+l) 

E  =  H  (B  +  —  +  B  )  +  4  d 

o  o  h..  o  2  o 

N 


The  constants  A  and  C  define  rigid  body  translation  of  the 
o  o 

laminate  as  a  unit.  The  remaining  constants  in  (35)  can 
all  be  expressed  in  terms  of  and  C^.  Hence,  two  constants 
which  effect  the  stress  distribution  have  been  introduced  in 
the  repeated  zero  part  of  the  homogeneous  solution. 

The  remaining  portion  of  the  complementary  solution  consists  o 
functions  of  the  form  (34)  corresponding  to  the  12 (N+l) -2 
nonzero  values  of  X  (we  are  assuming  that  these  roots  are 
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all  distinct) .  In  the  present  formulation  the  number  of  nonzero 
values  of  X  is  higher  than  that  obtained  in  [1] .  This  is  due  to  the 
difference  in  equilibrium  equations  for  the  global  domain.  The  requi¬ 
site  number  of  extra  boundary  conditions  are  obtained  from  the  last 
two  terms  of  (15) .  These  roots  occur  in  pairs  of  complex  conjugates 
a  +  ib.  Using  the  eigenvalues  and  eigenvectors  and  the  edge  condi¬ 
tions  (29  and  30)  we  can  obtain  the  solution  for  the  13(N+1)  func¬ 
tions  appearing  in  (19,  20,  23  and  25)  as 


12 (N+l)-4 


^  p(k)  eXmx  +  f(k)  +  f(k)f  (Jc=l , 2  ,N+1 ) 


m  =  1 


Ho  p 


where  the  last  two  terms  are  defined  by  (31-33,  35  and  36). 


The  force  and  moment  resultants  can  be  computed  by  substituting 
the  results  of  (37)  into  (17)  and  (18)  and  thence  into  the 


constitutive  relations. 


When  each  layer  is  isotropic  and/or  oriented  at  an  angle  of 
0°  or  90°,  the  compliance  components  S,c,  S_. ,  S and  S._ 

lb  Zb  Jb  4d 

and  expansional  strain  e^  vanish  in  each  layer.  This  leads  to 
the  vanishing  of  U(k),  <j>(k),  t^k)  ,  ,  N^k)  and  .  Conse¬ 

quently,  the  number  of  field  equations  and  boundary  conditions  is 
reduced.  This  case  must  be  treated  separately  by  specializing 
the  present  derivation  as  in  [1]. 

In  the  model  presented  in  [1],  laminates  consisting  of  a 


moderate  number  of  layers  (N>6)  could  not  be  analyzed  because 
of  computer  over flow/ under flow  limitations.  Furthermore,  while 
the  present  work  has  focused  on  a  fairly  special  case  (only  one 


global-local  interface) ,  more  general  arrangements  of  global  and 
local  domains  can  be  treated  by  simple  modifications  of  the  gen¬ 
eral  relations  given  here,  e.g.,  equations  (13a-13c)  must  reflect 
the  proper  positions  of  global  and  local  media.  In  this  way 
arbitrary  layers  can  be  modeled  in  a  local  fashion  to  define  the 
stress  field  in  the  entire  body,  if  desired. 
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SECTION  VI 

RESULTS  AND  DISCUSSION 

For  the  computation  of  numerical  results,  T300/5208 
graphite  epoxy  material,  with  the  following  elastic 
properties,  has  been  considered: 

Eu  =  20  X  106  psi,  E22  =  E33  =  1.4  X  106  psi 

G^2  =  G^3  *  0.8x10^  psi,  G2,  =  0.6  x  10^  psi 


V12  =  V13  =  °‘3  ' 


v 


23 


0.6 


where  E,  G  and  v  stand  for  Young's  modulus,  shear  modulus 
and  Poisson's  ratio,  respectively.  The  subscripts  denote 
the  corresponding  directions,  where  1,2,3  stand  for  x,y,z, 
respectively,  and  is  the  Poisson  ratio  measuring  strain 
in  the  j  direction  caused  by  uniaxial  stress  ck  . 

In  most  of  the  earlier  investigations  on  edge  effects, 
the  Poisson  ratios  v.^,  V13  anc^  v23  were  ta^en  to  be  equal. 

A  recent  experimental  study  has  revealed  the  values  shown 
above.  In  particular,  the  magnitude  of  v23  was  found  to  be 
equal  to  approximately  0.6  [16]. 

Figures  3-10  depict  the  distribution  of  stress  components 
o  ,  t  and  t  along  the  width  for  various  laminates.  The 

2  xz  y  Z 

abscissa  is  the  laminate  width  coordinate  normalized  by  the  half 
laminate  thickness,  such  that  in  these  diagrams  X=1  represents 
the  free  edge  of  the  laminate  and  X=0  represents  a  point  at  a 
distance  equal  to  the  half  laminate  thickness  from  the  edge. 
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I  1 

width  approaches  infinity  and  can  be  shown  to  be  very  accurate 

|  for  laminates  in  which  the  width  is  more  than  (approximately) 

\ 

;  twice  the  total  thickness.  The  coordinate  axes,  stacking  sequence 

and  loading  conditions  are  shown  in  the  figures.  In  the  symbolic 
notation  for  laminate  ply  orientations,  a  numeral  followed  by 
H,  Q  or  T  denote  one  half,  one  quarter  or  one  third,  respectively, 
thickness  of  the  corresponding  layer.  The  layers  under  a  bar 
constitute  the  global  region. 

Figure  3  shows  the  stress  component  oz  at  the  laminate  mid 
surface  versus  X  for  a  (0/+60)g  -  laminate,  calculated  by  using 
the  formulation  of  [1],  for  three  different  layer  thickness 

representations.  The  first  representation  is  such  that  the  I 

first  layer  from  the  midsurface  (0°)  is  modeled  as  two  sublayers, 
each  of  half  the  layer  thickness  h/2,  and  the  other  two  layers 
are  treated  individually.  In  the  second  representation,  the  0° 
layer  is  modeled  as  three  different  sublayers  of  equal  thickness 
h/3.  The  third  representation  considers  the  0°  layer  as  three 
sublayers  of  h/3  thickness  each  and  the  60°  layer  as  two  sub¬ 
layers,  each  of  thickness  h/2.  The  fact  that  the  results  by  all 
the  three  representations  of  the  laminate  are  the  same  show  that 
■  these  results  are  nearly  exact. 


Figure  4  shows  the  stress  component  o for  (0/(+60)2)s  - 
laminate  by  three  different  representations.  For  the  first 
representation,  the  theory  developed  in  [1]  has  been  used  and  for 
the  other  two  representations  the  present  global- local  model  has 
been  used.  It  has  been  seen  that  the  results  obtained  by  the 
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middle  representation  i.e.  (OH  /OH  /60/-60/60/-60)  are  nearly 
identical  with  those  obtained  by  [1]  for  the  same  representation 
of  the  local  domain.  The  values  obtained  through  the  third 
representation  only  differ  slightly  from  those  for  the  other  two 
cases.  Further  analysis  shows  that  the  "hump"  is  caused  by 
insufficient  subdivision  of  0°  layer,  rather  than  the  local- 
global  model. 

Figure  5  shows  the  variation  of  the  stress  component  c?z 
along  x-axis  for  (90H  /+30/90 ) s~laminate  by  three  different 
representations.  It  is  seen  that  the  results  obtained  by  the 
local-global  models  differ  considerably  from  those  obtained  by 
model  of  [1].  The  results  by  the  third  representation  differ 
from  those  for  the  second  representation  adjacent  to  the  free 
edge  and  are  close  elsewhere.  Figure  6  shows  more  results  for 
the  same  laminate.  In  this  figure  the  third  representation  is 
the  same  as  that  of  Figure  5.  Also  in  Figure  6  we  have  used  a 
different  representation,  (90Q/90Q/-30/30/90 ) ,  of  the  same  laminate. 
In  this  case  it  has  been  found  that  the  agreement  between  the 
[1]  model  results  and  the  local-global  model  results  is  again 
quite  good.  Hence,  an  extra  hump  in  the  results  by  the  third  rep¬ 
resentation  is  likely  due  to  insufficient  subdivision  of  the  inner 
90°  layer.  This  shows  that  another  factor  may  be  important  in 
obtaining  satisfactory  results.  The  precise  definition  of  this 
factor  is  not  known  presently,  however,  it  appears  that  a  gradual 
transition  between  the  local  and  global  regions  may  be  helpful  in 
obtaining  accurate  results,  i.e.  the  middle  representation  of 
Figure  6.  Another  factor  involved  in  the  results  of  Figure  5  is 
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3  computed  stress  component  a  is  of  considerably  lower 
de  as  compared  to  that  in  the  laminate  of  Figure  4.  Thus, 
olute  magnitude  of  the  error  in  the  results  for  the  first 
ntation  of  Figure  5  is  small,  although  the  relative  error 
quite  large.  It  can  be  seen  from  Figure  7  that  for  the 
.y  representations  as  those  in  Figure  5  the  relative  error 
results  for  (90H  /+15/90)  -  laminate  is  small  as  compared 

:  for  (9 OH  /+30/90)g  -  laminate.  The  results  for  (90Q/ 
5/15/90)  representation  are  also  computed  and  are  the 
s  those  of  the  first  representation  of  Figure  7.  Thus, 
ears  that  the  absolute  magnitude  of  the  stress  component  cz 
o  important  in  modeling  the  laminate  representations  for 
te  results. 

'igure  8  shows  the  stress  component  txz  at  the  90/- 30  interface 

!90H  /+30/90 )  -  laminate.  As  in  the  case  of  a  ,  the  results 

s  z 

lei  [1]  are  nearly  identical  to  those  by  the  present  model 
[90Q/90Q/- 30/ 30/90)  representation.  The  values  computed 
(h  the  other  representation  of  the  present  model  are 
:ly  different  from  the  others.  However,  the  maximum  value 
will  increase  with  larger  numbers  of  sublayers  in  the  90° 

This  is  a  characteristic  of  the  general  class  of  models 
presented  and  is  discussed  in  more  detail  in  [1]. 

’igure  9  gives  the  stress  component  x^z  for  the  afore- 

>ned  laminate  at  the  90/-30  interface.  In  this  case  also,  the 

■ison  amongst  the  results  by  three  different  representations 

tsonable,  although  an  elastic  singularity  is  expected  in 

(tress  component.  Hence,  again  significant  dependence  on 

rer  size  will  be  present  near  the  edge. 
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Figure  10  shows  the  stress  component  o for  (0/(±60)n) 
n=l,  2,  4  laminates  as  computed  by  the  lo~al-global  model.  The 
results  given  for  the  first  two  values  of  n,  i.e.,  n=l,  2,  are 
already  shown  to  be  identical  with  those  obtained  by  using  [1] , 
Figures  3-4.  The  results  for  n=4  show  the  expected  trend.  There 
exist  no  results  in  the  open  literature  to  compare  with  these 
values . 
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SECTION  VII 
CONCLUSIONS 

A  self-consistent  model  has  been  developed  to  investigate  the 
stress  fields  in  laminated  media  consisting  of  numerous  layers. 

The  new  model  defines  detailed  response  functions,  such  as  inter¬ 
laminar  stresses  and  single  layer  forces  and  moments  in  a  pre¬ 
determined  region  of  interest  (local) ,  while  the  remainder  of  the 
domain  is  represented  by  its  effective  material  properties  and  the 
corresponding  resultant  forces  and  moments  (global).  The  local 
model  employs  a  theory  [1]  which  approaches  the  theory  of  elasticity 
in  the  limit  of  vanishing  layer  thickness.  The  global  model  is 
based  upon  the  theory  given  by  Whitney  and  Sun  [6]  which  has  been 
demonstrated  to  produce  good  agreement  with  elasticity  results  on 
the  global  boundary  for  a  particular  laminate  by  Pagano  [7].  While 
a  particular  arrangement  of  global  and  local  domains  has  been 
considered  here  for  brevity,  there  is  no  difficulty  in  extending 
these  results  to  include  more  general  arrangements,  including 
the  use  of  more  than  one  global  domain.  The  importance  of  the 
latter  option  follows  from  the  observation  that  model  accuracy 
may  be  improved  by  a  gradual  rather  than  abrupt  transition  of 
region. 

The  effectiveness  of  the  model  has  been  demonstrated  by  use 
of  numerical  examples  based  upon  the  free-edge  class  of  boundary 
value  problems  in  laminate  elasticity.  Preliminary  results  have 
been  shown  to  be  very  promising  although  an  apparent  loss  in  accuracy 
occurs  in  the  calculation  of  stress  components  of  small  magnitude, 
which  may  thus  require  finer  subdivision  of  the  local  region  than 
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would  normally  be  required.  Similarly,  the  effect  of  the  afore¬ 
mentioned  "transition  region"  on  model  precision  will  require 
further  study.  These  studies  as  well  as  the  development  of  a  solutio 
schemes  for  fully  three  dimensional  problems  (which  will  depend 
only  on  two  space  variables  in  this  theory) ,  will  be  subject  of 
future  investigations. 

It  is  clear  that  theories  of  the  type  presented  here  are 
needed  to  describe  the  response  of  laminated  structural  components 
used  in  practice.  However,  experimental  activity  in  this  regard 
is  vital,  as  proper  interpretation  of  the  field  analysis,  parti¬ 
cularly  in  regions  of  very  steep  stress  gradients,  is  needed  to 
characterize  initial  failure  and  subsequent  damage  growth  in  these 
bodies . 
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Figure  la.  Laminate  Half  Thickness  Divided  into  more  than  One 
Global  Domain  with  Different  Types  of  Interfaces. 

I  -  Local-Global  Interface,  II  -  Global-Local  Inter 
face.  III  -  Global-Global  Interface. 


Figure  lb.  Laminate  Half  Thickness  with  One  Local-Global 
Interface . 
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3.  Stress  Distribution  c2/(exl0  )  PSI  versus  Width 
Coordinate  X  at  the  Mid  Surface  of  the  Laminate 
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Figure  4.  Stress  Distribution  o^/(ex 10  )  PSI  versus  Width 
Coordinate  X  at  the  Mid  Surface  of  the  Laminate. 
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Figure  5.  Stress  Distribution  o?/(exlOD)  PSI  versus  Width 
Coordinate  X  at  the  Mid  Surface  of  the  Laminate. 
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Figure  6.  Stress  Distribution  oz/(exlO  )  PSI  versus  Width 
Coordinate  X  at  the  Mid  Surface  of  the  Laminate 
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7.  Stress  Distribution  oz/(cx 10  )  PSI  versus  Width 
Coordinate  X  at  the  Mid  Surface  of  the  Laminate 
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figure  8.  Stress  Distribution  tx-/(ex10^)  PSI  versus  Widt 
Coordinate  X  at  the  90/-30  Interface  of  the 
Laminate. 
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Figure  10.  Stress  Distribution  o^/(exl0  )  PSI  versus  Width 
Coordinate  X  at  the  Mid  Surface  of  the  Laminate. 
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